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Abstract 

We consider magnetic catalysis in a field-theoretic system of (3+l)-dimensional 

Dirac fermions with anisotropic kinetic term. By placing the system in a strong ex- 

_ , ternal magnetic field, we examine magnetically-induced fermion mass generation. 

D I When the coupling anisotropy is strong, in which case the fermions effectively lo- 

'^ ' calize on the plane, we find a significant enhancement of the induced mass gap 

compared to the isotropic four-dimensional case of quantum electrodynamics. As 

^ ■ expected on purely dimensional grounds, the mass and critical temperature scale 

H , with the square root of the magnetic field. This phenomenon might be related to 

recent experimental findings on magnetically-induced gaps at the nodes of d-wave 

superconducting gaps in high-temperature cuprates. 






1 Introduction 

The phenomenon of magnetic catalysis, first suggested in [|l|, and subse- 
quently developed also in 0, namely the dynamical generation of a fermion 
mass gap in the presence of external magnetic fields, has wide applications, 
ranging from particle physics and physics of the Early Universe to con- 
densed matter Q, |^, |^ . 

In the case of four- dimensional Abelian f/(l) gauge- Dirac-fermion mod- 
els the magnetic catalysis phenomenon, for an external magnetic field B, is 
known to yield a relatively small mass gap m^d = m^diB) (for standard 
low-energy electromagnetism, with fine structure constant of order 1/137) 
and an associated critical temperature of the same order Tg ^ rn^d- Such 
analyses pertain to Abelian gauge interactions with spatially isotropic cou- 
plings. It is the purpose of the present article to consider spatially anisotropic 
gauge interaction, in such a way that the gauge coupling on a spatial plane, 
defined by, say, the x, y directions of a three-dimensional space, is much 
stronger than the coupling along the ^-direction. 

Physical motivation for such studies is provided by recent experimen- 
tal findings in the physics of high-temperature layered superconductors [§]. 
According to such experiments, one looks at the thermal conductivity proper- 
ties of a sample of high-temperature superconductors in the superconducting 
phase, in the presence of strong external magnetic fields (with intensities up 
to a few Tesla). The high-temperature superconducting samples are known 
to be strongly type II, and so the magnetic field lines penetrate the material 
significantly. Moreover, such materials are known to be d-wave superconduc- 
tors with a layered structure of Cu-0 planes, characterized by nodes in their 
superconducting excitation gaps. The experiments of ||^ have demonstrated 
that below a given temperature, which is smaller than the critical tempera- 
ture of the superconductor, and which scales as the square root of the applied 
magnetic field, there are plateaux in the thermal conductivity diagrams. The 
plateaux indicate the opening of gaps at the nodes below this new 'critical' 
temperature (not to be confused with the superconducting/normal phase 
critical temperature) . 

In ref. Q] we suggested that this phenomenon might be a straightfor- 
ward application of the magnetic catalysis phenomenon |l|, |], ^ to the three- 
dimensional effective problem of physics near the nodes. Indeed, lineariz- 
ing the excitations about such nodes one obtains, upon the assumption of 



spin-charge separation [|iy], a relativistic electrically-charged ferniion (Dirac) 
system, coupled to statistical gauge fields in the presence of external elec- 
tromagnetic fields. The strong external magnetic field induces the opening 
of a nodal holon gap, which scales with the magnetic field B. The exper- 
iments of PI have shown a square root scaling with B. In the models of 



[^, 0, the statistical gauge fields represent effective spin-spin antiferromag- 
netic interactions, which are believed by many to be relevant for the physics 
of the high-temperature superconductivity. In ref. 0], and in all the subse- 
quent works [^, ^, where various interactions, including four-fermi, among 
the holons have been considered, the model systems have been assumed to live 
exclusively in (2 + l)-dimensions, ignoring completely any four- dimensional 
physics effects. 

This may not be physically correct, especially from the point of view 
of the gauge and electromagnetic interactions, which are known to be fully 
four- dimensional Q. Moreover, the superconducting gaps are actually four- 
dimensional, although strongly anisotropic, i.e. with dominant components 
along the Cu-0 planes, but with suppressed, however non- vanishing, compo- 
nents along the perpendicular (interplanar) direction. 

It is the purpose of this article therefore to consider such an anisotropic 
four- dimensional situation and study the consequences for magnetically in- 
duced dynamical mass generation for the fermionic excitations. From the 
condensed-matter point of view, such excitations may be the (continuum 
limit of) holons, carrying the electric charge only, but no actual spin (in 
the microscopic sense). In ref. [^] the holons were assumed purely three- 
dimensional (planar), as a result of the localization of their wavefunctions on 
the superconducting planes. This is a basic feature assumed to characterize 



the microscopic physics behind spin-charge separation ||T^. In the present 
article the fermions will be assumed four-dimensional, but with anisotropic 
kinetic terms, which allows for interplanar hoping. In the strong anisotropic 
case one should recove the three-dimensional case. In this sense one should 
have the extension of the concept of spin-charge separation to four dimen- 
sions, in the strongly anisotropic case. This, of course, does not apply to 



the concept of fractional statistics |11], which is only an exclusive feature of 



"'^For instance, interlayer couplings via magnetic spin-spin interactions are known to 
exist in the planar high-temperature superconducting materials, and, of course, the elec- 
tromagnetic interactions are fully four dimensional. 



the three-dimensional case. However, there may be a different scenario, in 
which the fermions discussed here are viewed as related to real nodal electron 
excitations (carrying both spin and charge), capable of interlayer hopping. 
It is because of this latter interpretation that, when we discuss in section 
5 the non-relativistic situation of relevance to condensed matter, we shall 
ignore any statistical gauge interactions among the fermions, keeping only 
the electromagnetic interactions. 

In case one assumes the existence of four-dimensional holons, by ex- 
tending the concept of spin-charge separation to four-dimensions, allowing, 
though, for weak interlayer hopping of holons, one might encounter a sit- 
uation relevant to the fully relativistic case discussed in the next section, 
in which the gauge fluctuations represent statistical gauge fields. In such 
a case one may even encounter non-abelian gauge interactions. Indeed, in 
the physical models of []TT|, the statistical gauge fields are non-Abelian, of 



SU{2) type, due to an underlying 'particle-hole symmetric' formulation of the 
spin-charge separation ansatze. Such interactions are not responsible for the 
opening of holon mass gaps, because of the fermion spectrum in the three- 



dimensional models of |jTl[. It is the statistics-changing ?75'(l)-interaction, 
exclusive for three-dimensional systems, which is strong enough to generate 
such gaps. Unfortunately the Us{l) interaction does not have an analogue 
in four- dimensions. Thus, although probably relevant for the superconduc- 
tivity scenaria, nevertheless such statistical interactions may be not directly 
relevant for the four-dimensional physics underlying the findings of [§]. 

For our purposes below, therefore, we shall ignore such non-Abelian sta- 
tistical interactions, and concentrate rather exclusively on the role of real 
electromagnetic interactions, described by a potential A^, in inducing, under 
the influence of a strong external magnetic field, a mass gap for the nodal 
holons. We shall assume, however, that the electromagnetic interactions are 
screened along the z-direction (interplanar) , and in this sense we shall treat 
the gauge coupling e as spatially anisotropic. A subtlety of this model is the 
relativistic non-invariance of the electromagnetic interactions in the presence 
of holons. The fermion (holon or, even, charged electron excitations) part is 
a relativistic system describing nodal excitations, but for such a system the 
role of the limiting speed of 'light' is played by the fermi velocity vp of the 
nodes. On the other hand, the real electromagnetic field propagates with the 
velocity of light c':$> vp (for realistic systems c ~ lO^f f, but for us vp will be 
considered as a phenomenological parameter). When the combined system 



is expressed in terms of the fernii velocity, the electromagnetic interaction 
will be non relativistic, with the velocity of light c appearing explicitly in the 
Maxwell action. This leads to a non-relativistic form for the photon propa- 
gator in the environment of holons or equivalently to a non-relativistic form 
for the fermion propagator if we rescale the fields. 

Our approach in the present article will be field theoretic, and we shall not 
attempt to make further contact with the condensed-matter systems, apart 
from the very generic features mentioned above. However, as we shall see, 
our findings are interesting enough, and indeed may be of use in attempts to 
explain the phenomenon of [§] by means of the magnetic catalysis mechanism 
of relativistic fermions [§, ^, ^. A brief discussion on such 'phenomenology' 
will appear at the end of our article. 

At this stage the reader's attention is drawn to a very interesting recent 
work ll^ , which analyses the phenomenon of chiral symmetry breaking on a 



brane domain wall, embedded in a higher- dimensional space time, in the ab- 
sence of any external fields. This situation is different from the one discussed 
here, not only because there is no external field, but also because in the work 



of [0 the fermions are completely localized on the brane, while the gauge 
fields propagate in the bulk. In contrast, in our case we allow interlayer hop- 
ping for fermions, which is anisotropic, and also we have anisotropic gauge 



couplings. In the scenaria of |T2[, the fermions could represent purely three- 
dimensional holons of planar high-temperature cuprates, localised on the 
Cu-0 planes. In such scenaria, therefore, the spin-charge separation would 
be an exclusive feature of the planes, not extended (even in the strongly 
anisotropic sense) to four dimensions, but the statistical interactions (repre- 
sented by the gauge fields) could extend to interplanar coupling situations. 
This would make an important physical difference from the scenaria discussed 
in the present article. 

The structure of the article is as follows: in section 2 we study the mag- 
netic catalysis phenomenon in relativistic gauge-fermions systems in the pres- 
ence of anisotropic four- dimensional gauge couplings. In section 3 we demon- 
strate that, in the strongly anisotropic case, there is a significant enhance- 
ment of the induced mass gap nidyn (on the plane), as compared with the 
isotropic case. In section 4 we study the system at finite temperature, and 
compute the critical temperature T^ above which the magnetically-induced 
mass gap disappears. We demonstrate that Tc ~ m^y„, as expected on natu- 
ral grounds. The square-root scaling with the magnetic field intensity is also 



demonstrated. In section 5 we attempt to make contact with realistic con- 
densed matter systems, and the experiments of [§, by repeating the above 
analysis but for a non-relativistic fermion system, coupled to a relativistic 
electromagnetic field, again with anisotropic couplings, whose quantum fluc- 
tuations are taken into account in the absence of any other interactions. 
The role of the limiting speed of light for the fermions (holons or electrically- 
charged nodal excitations) is played by the fermi velocity vp, while, the speed 
of light is c > vp- In this case the ratio vp/c < 1 is a phenomenological pa- 
rameter of the model. The main features of the previous fully-relativistic 
case are maintained, especially as far as the enhancement of the fermion gap 
is concerned in the strongly anisotropic case. However, there are extra sup- 
pression factors by powers of vp/c, relative to the relativistic case of section 
2. Conclusions and some discussion, with relevance to 'phenomenology' of 
high-temperature superconductivity, are presented in section 6. Some formal 
aspects of the Schwinger- Dyson analysis for the anisotropic case are presented 
in an Appendix. 

2 Anisotropic gap equation 

The Lagrangian density which includes the anisotropy is, in the absence of 
an external field 

C = -\f^,F^^ + ^ [, ^ _ e yl - X (2^3 - eAg) 7' - m] V-, (1) 

where ^ = 7^9^, /i, z/ = 0,1,..., 3, 7^ are four-dimensional 4x4 Dirac 
matrices, ip are four- component spinors, and the parameter x controls the 
anisotropy. The case x = corresponds to a totally isotropic situation (the 
usual quantum electrodynamics (QED)), while x = 1 corresponds to a to- 
tally anisotropic one: the gauge field lives in 3-1-1 dimensions whereas the 
fermions are effectively localized in 2-1-1 dimensions. We choose the 7 matri- 
ces 7*,z = 0, 1,2 diagonal and the matrix 7^ non-diagonal, as in 0. In the 
totally anisotropic case a; = 1, 7^ does not appear in the action and the other 
7 matrices decompose in two sets of two-component Dirac matrices, appro- 
priate for the irreducible Dirac algebra in (2 -|- l)-dimensions. In such a case 
the four- component spinors also decompose to an even number (2 for a single 
flavour) of two-component (2-|-l)-dimensional spinors, and one recovers the 



planar case, after integration over the third dimension |. 

It should be stressed that Eq.(P of course respects the gauge invariance. 
With such a Lagrangian, and in the absence of external field, the bare fermion 
propagator is given by 

iS^'^{p) = ij) — xp'i'y^ — m. (2) 

As an illustration of the effect of the anisotropy, we can compute the propa- 
gation of a fermion in the direction 3. Let us define in the Euclidean space 



(j){r-i) = tr dr4dr2dr I < 0\ilj{0)ilj{r)\0 > 

= tr/^5(0,0,p3,0)e^^^^^ (3) 

J ZTT 

i.e. the propagation rate of the fermions from the 2+1 dimensional worlds 
r3 = and r^ ^ 0. A straightforward computation leads to 

^(,,) = _exp^-^j. (4) 

which shows that as the anisotropy increases, the fermion propagation in the 
direction 3 decreases exponentially. Eventually, as x ^ 1, this propagation 
vanishes if r^ ^ 0. This corresponds to an effective delta-function in r^, which 
kills the r^ integration in the action leading to a three-dimensional theory. 

We derive in the appendix the Schwinger- Dyson equation for the fermion 
propagator, taking into account the anisotropy parameter x and find (we do 
not write the space-time indices) 

G = S-47Ta f S-f^'GA'^GD^^ + xAna f S-f'^GA^GDs^ (5) 

where G is the full fermion propagator, D^,^ the full photon propagator and 
A'' the full vertex which satisfies at the tree level 



^In the physical case of relevance to condensed matter, the third dimension may extend 
between two Cu-0 layers, which play the role of boundaries of the available space. In such 
a case the extra (third) dimension is integrated in the bulk space between these boundaries. 
In the analysis below we shall not assume explicitly such geometries, but rather treat the 
anisotropic four-dimensional case in a generic sense, specified by the Lagrangian p). 
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K,, = r if /i ^ 3 

AL. = (l-x)y (6) 

We use the usual definition of the fine structure constant: e^ = Ana. 

In the presence of an external, constant and homogenous magnetic field 
in the direction 3, we can choose an external gauge such that Af^^ = for 
which the parameter x will not couple to the external field and thus will only 
play a role in the substitution p^ -^ {l — x)ps. The lowest Landau level (LLL) 
approximation [^j for the fermion propagator will be then 

S^{y, z) = e^<=?''^^rW5i(y _ z), (7) 

where the Fourier transform of the translational invariant propagator S^ is 

S^p) = e-'i'l"" ,'^J-"'y\ (8) 

where p^ = p\+p\ is the transverse momentum squared. We will take m = Q 
for the bare propagator (not taking into account the interaction with the 
dynamical gauge field) and m = m{po,p3), the dynamical self energy which 
depends only on the longitudinal momenta in the LLL approximation, for 
the full propagator. The integral equation describing the magnetic catalysis 
in the LLL approximation contains only the components Dqq and D33 of 
the photon propagator [^, due to the spin projector (1 — i7^7^)/2 in the 
fermion propagator (Q). These two components lead to equal contributions 
and thus the second integral in the Schwinger-Dyson equation (^ which is 
proportional to x gives half the contribution of the first integral. The integral 
equation reads then, if we neglect the corrections to the vertex. 



G^ = S^ -47Tall-x + —] f 5S"G'S"G'^^|j,|| 



(9) 



where || denotes the longitudinal components. The final integral equation is 
found by making the substitution ps — i> (1 — x)p5 in the fermion propagators 
appearing in the integral of Eq. (^ which has already been computed in the 
isotropic case ^. If we rescale all the quantities homogenous to a mass by 
eB\, we have in the Euclidean space 
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fJ'ih, h) = \l-x + 



2 / 7r2 7 (^i-xypl+pl + fl'^{p3,p4) 

""Jo 2u+{p,-hf + {p,-hr ^^°^ 

where /i is the dhiiensionless self-energy and u a dimensionless transverse 
momentum squared. The momenta ps and p^ play a symmetric role in 
the photon propagator whereas they enter non-symmetrically in the fermion 
propagator, which leads to a quite difficult analysis of the integral equation 
if we wish to take into account the momentum dependence of the fermion 
self energy. We know that this momentum dependence is essential in 3+1 
dimensions [§ but not in 2+1 dimensions, as long as we consider a qualita- 
tive description [|I^. In this paper we wish to give a qualitative description 
of the magnetic catalysis in the strongly anisotropic regime 1 — x << 1 and 
thus will make the constant self-energy approximation in which the integral 
equation finally reads 



x^\ a f dp3dp4 r°° due " 

2 / 7r2 y (1 - xfpl + pI + nl h 2u + pI+pI' 



or, when we make the angular integration. 



x^\ a f°^ dp r°° due " , , 

l-a; + — -/ ^ : / 12 

^J^Jo ^(p + ;,2)p_^)2p + ^2]yo 2u + p 



The analysis of the approximation ([Tl|) in the isotropic regime x << 1 will 
not lead to a reliable quantitative dynamical mass /xq but will show us the 
qualitative tendency of /xq to increase when the anisotropy increases. 

3 Enhancement of the dynamical mass 

In this section we will show that the anisotropy generates a considerably 
enhanced mass gap, compared to the isotropic situation where we know that 
in the constant self-energy approximation [|1|], P], 



/io(x = 0)~V2exp -,/- . (13) 

V V a/ 

To see the increase of the dynamical mass when x increases, we make an 

expansion of the integral equation ([TTD up to the order x^ and find 

1 = ^ ^" due-"" |ji + X (2J2 - Ji) + x^ (^4X3 - 3X2 + ^Xi) +...}, (14) 



where 



Xi 



dpsdpi 



TC 






X. 



2u — /io 

dpsdpipl 



X, 



(P3+P4 + /^o)'(P3+P4 + 2m) 

TT 2m ln(2M//io) — 2m + /^o 
2 (2« - /ig)2 

dpsdpApt 

(P3+P4 + /io)'(P3+P4 + 2M) 

37r 8^2 \n{2u/nl) + (/"o " 6m) (2m - /ig) 
16 (2m - /ig)3 ■ 



(15) 



Despite their appearance, the integrals ([15D are converging when 2m — >■ /Zq, 
what can be checked by expanding the logarithms around 2u = /ig- 

For the integration over the transverse momentum u, we write that for 
/io << 1: 



., ln(2M//ig) 

2m — yUo 

TT /"^/aiq , Inw vr /"^/aiq _ Inw 



/•CXD /"OO 

/ dulic^'^ = TT / (iwe 

JO JO 



^*_^ 



(iw ~ — / (iw 



2 JO M — 1 2 Ji u 



2/'v^' 



vrlnM^ . (16) 

V A'o / 



Similar approximations lead to 



2 \ /io y 



/o 8 \ /io 

and give the following equation for the dynamical mass 



r,,a,e- . |ta^(^). (17) 

Jo 8 \ /in / 



(-^■■■)-(S)4- 

We have then for the dimensionful dynamical mass, to order x^, 

m,,„(x)~v^V^exp|-y^h -^U. (19) 

We see then that the anisotropy {x > 0) has the effect to increase the dy- 
namical mass. 

Let us now come to the totally anisotropic regime and compute the dy- 
namical mass obtained for x = 1. In this case, we find from (|1^) 



l = ^/_ ^SL=( ,,i!^^. (20) 



a r°° dp r°° due " 
27r io Vp + 1 io 2m + ppq 

To perform the integration over p, we suppose that 2m > Pq, which consists 
in neglecting the interval [0, Po/2] in the integration over w, which will be 
justified since we will find Pq << 1. We have then 

1 ~ / rfwe'" , = - - / du , (21) 

2/^0 ^g/2 J2u - pg 4 7i Vm-1 



such that we finally obtain for the dimensionful dynamical mass 



m^ynix = 1) ~ —=s^\eB\. (22) 

This result shows that the anisotropy plays a fundamental role in the gen- 
eration of a big mass gap since ra^ynix = 1) >> mdyn{x = 0). We show 
in figure |I| the dynamical mass as a function of x, for a given value of the 
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Figure 1: /xq versus x for a = .01 



coupling a. This figure has been done with the numerical study of Eg . (|T^) 
and we note the logarithmic scale, showing the exponential increase of the 
dynamical mass. The analytical results (0) and ( p2D are confirmed by the 
curve in the asymptotic limits x << 1 and 1 — x << 1. 



4 Critical temperature in the anisotropic regime 

We now turn to the finite temperature treatment of the problem so as to find 
the critical temperature when x = 1. 

The integral equation at finite temperature is obtained by the usual sub- 
stitutions "p^ ^* ijJi = {21 + l)7rt and J dp^ -^ 27itJ2i ^^ the equation ([TT|), 
and t = T / J\eB\ is the dimensionless temperature. We just have to pay 
attention to the photonic momentum p^ which has to be changed by uji—ujo, 
that is to say to a bosonic Matsubara mode. We obtain then the following 
equation for the critical temperature defined by ^o{tc) = 0: 
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o- 


/■ , 1'°° 




dp J 


TT J 


Jo 



l=oo -1 

due-% y ^n7^ 5 7 ^ (23) 

where ujq = vrtc- The summation over the dimensionless Matsubara modes 
uji is done with the usual contour deformation trick detailed for example in 
[|l^ and gives 



'^ 1 11 



+ , ' coth^V^" + "^^ -l-^u-pl 



2J2u + pl \ 2tc ) {ujI + 2u + pI) 



,2^2 • 



We suppose that the coth is close to 1 since tc << 1- This is actually not valid 
when 2u + p1 « 1 but this region is negligible in the remaining integrals. 
The integral equation reads then 



^r^^e-o/^- ^ 



^Jl + u 

+ ^ r due-y^ r dv ^-"-"; (25) 

7c Jo Jo {1 + u + v^yVu + v^ 



The dominant contribution in Eq. ( PSJ ) comes from the first term when t^ « 
1 since in the second term, the integration over v gives 



°° I — u — v^ 

dv- 



{l + u + v'^y^/u + v"^ 

1 , 1 , /2 + n + 2VTT^\ ^^^ 



l + u 4(l + n)3/2 ^2 + u - 2y/TTu J ' 
which leads to higher powers of tc- The critical temperature is then given by 

1 ~ — / du , (27) 

such that the dimensionful critical temperature is finally 
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Tc{x = 1) ^ ^V^ = l^ciynix = 1). (28) 

This result confirms what was expected, i.e. the critical temperature is of 
the order of the dynamical mass at zero temperature. The same conclusion 
was found in M for isotropic QED. 



5 Non-Relativistic Gauge Field— Fermion 

Models and an Application to Condensed 
Matter 

So far we have examined a relativistic system, in which the gauge interactions 
may be assumed distinct from the real electromagnetic interactions. From 
a condensed matter view point, such systems may have some relevance to 
effective gauge theories of t — j models, with intrasublattice hopping, per- 
taining to spin-charge separation scenaria applied, however, to nodal (i-wave 



excitations 14 



In this section we shall consider a problem which might be considered 
as more "realistic" , in the sense of being connected directly to the observed 
physics of high-temperature superconductors, in view of the recent experi- 
mental findings of . We shall not assume ad hoc any statistical interaction, 
but we shall continue to apply a spin-charge separation scenario |jlO| for the 
nodal excitations of the high-temperature (i-wave superconductors. 

In this case, the electrically-charged excitations (holons) around the nodes 
of a d-wave supercondcuting gap will be represented by (3 + 1) -dimensional 
Dirac fermions coupled to an electromagnetic field, described by a potential 
A^ = {AQ,Ai). The four-dimensional nature of the fermions implies the 
possibility of interplanar hopping for such excitations in the materials, which 
is a realistic feature. We shall consider the case where the electromagnetic 
field has a background, corresponding to a constant magnetic field along the z 
direction, and four-dimensional quantum fiuctuations around it, and assume 
no further interactions among the holons. However, the electromagnetic 
coupling e will be assumed screened along the z-direction, i.e. the coupling 
is anisotropic. Such a screening may be provided by a combination of both 
the chemistry and geometry of the planar materials, which consist of layers 
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of Cu-0, with lattice spacing of 5 Angstroms, while the interlayer distance is 
much larger, of the order of 100 Angstroms. In the interlayer space there are 
doping atoms etc., which may be responsible for an effective charge screening. 

From a formal view point we might think of the effective electric-charge 
anisotropy as follows: the coupling constant e'^/Air entering the problem, as 
a result of quantum fluctuations of the electromagnetic fleld, is actually a 
running coupling constant, which depends on the available energy. If the 
latter is viewed as the inverse of a characteristic distance for the problem, we 
then observe that, as a result of the relatively large separation between two 
holons (which are the relevant degrees of freedom for the problem) at different 
layers, assumed perpendicular to the z direction, the resulting effective charge 
(along the z direction) lies in the infrared regime, and as such is much weaker 
than the effective charge between holons in the same layer and actually at 
characteristic distances of the order of the inverse of the magnetically- induced 
mass gap, which we are interested in for the purposes of the present work. 
In high-temperature materials this distance, i.e. the magnetic coherence 
length of nodal charged excitations, is short enough and of the order of a 
few angstroms, which justifles the assumed strong anisotropy of the electric 
charge. 

With these in mind, we consider the following model. The photon kinetic 
term is the usual one 

Cg = -\f^uF^\ (29) 

with do = d/dt (c = 1) and the free fermion kinetic term is 

CF = ^{t-f'^do + ti''dk)iJ (30) 

with Oq = d/d{vpt). Since the fermion density and current are (the notation 
A denotes spatial three vectors): 

p = tp^tjj 

i = VF^Pii^, (31) 

the interaction between the fermions and the gauge fleld is 

Ci = e^7°AoV^ + evF^l-M (32) 
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such that the final Lagrangian will contain the following fermionic part 



^ [r^ (^9o7° - ^Ao7°) + {^^kl^ - gAu^'') - x {id^,^^ - gA^^^)] ^, (33) 

where we define rj = l/vp, g = evp and k = 1,2,3. The bare fermion 
propagator is then, in the absence of external field, 

iS'\p) = rjpo-f^ + pkl^ - xp37^ - m (34) 

It is straightforward to see from the derivation shown in the appendix that 
the Schwinger-Dyson equation changes to 



G = S-r]47ra f S-f^GA^GDo^ 

-Ana f S-i^Gk^GOku + xAna f S-f^GA'^D^u, (35) 



where 



Ai;,, = Y if ^ ^ 0, 3 



tree 



ALe = (1-^)7', (36) 

and g"^ = Ana. The integral equation corresponding to (|TTp is finally 



h' + (1-^)^1^ ,, ■„?".„.,.. / „.. ::. . .r (37) 



^2-, a r dp^dpi r^ due ^ 

'2^7 (l-xypj + riYA + f^o-lo 2u + pI+pI 

This last equation reads for x = 1 

a f dp3dp4^ f°° due'"^ 



2tt^ J pI + (Ho/ijy Jo 2u + pI+pI 



^2' 



(38) 



which is the equation (|20| ) with the substitution /xq — > /io/'?- Since t] >> 1, 
the condition fio/rj << 1 is still valid and the result is then for the dimen- 
sionful dynamical mass 
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3/2 e 



2 



^dyn^Vjl ^^Vl^^l- (^9) 



For the finite temperature case, the equation ( ^31) does not change (the 



factors 7] cancel when x = 1) and the critical temperature is then 

o2 



e 



T^c^vT^-^\l\eB\ (40) 



It is important to notice that both the the dynamical mass (pQ]), and the 
associated critical temperature (^Oj), are proportional to a suppression factor 
Vp = O{10~^) and vj = O{10~^^) respectively, compared to the previoulsy 
considered relativistic case. This may have important phenomenological im- 
plications when one attempts to compare the scenaria advocated here with 
realistic condensed-matter situations ||^. We shall not do such analyses here, 
given that at present we lack a detailed derivation of the continuum models 
discussed above from microscopic condensed matter models. This is essential 
for providing the correct order of magnitude of the various coupling constants, 
such as gauge, hopping elements etc., entering the model. 

6 Discussion 

In this work we have discussed the case of four-dimensional anisotropic 
fermions coupled to anisotropic-coupling Abelian gauge fields, and exter- 
nal magnetic fields. By considering fluctuations of the gauge fields we have 
considered the magnetic catalysis phenomenon, i.e. the dynamical mass gen- 
eration for fermions under the influence of strong external magnetic fields. 
We have considered two cases: the fully relativistic case, in which fermions 
and fluctuating gauge fields are relativistic, and the case where the photon 
fields are relativistic, but the fermion part of the Lagrangian is non relativis- 
tic, with a fermi velocity vp ^ c {c is the speed of light) playing the role of 
the limiting velocity for the fermionic part. 

In both cases a strong planar anisotropy has been assumed. The analysis 
has shown a significant enhancement of the induced mass gap, relative to the 
isotropic four- dimensional case. In the non-relativistic case, however we find 
that the induced mass gap is found suppressed by some power of the fermi 
velocity vf < c, as compared to the relativistic case. In realistic situations, 
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vp/c ~ 10~^, and the suppression factors are significant. In both cases, the 
mass gap and the associated critical temperature, above which the magnetic 
catalysis disappears due to thermal disorder, exhibit a square-root scaling 
with the magnetic field intensity, as expected on natural grounds. 

The non relativistic case, may be related to the situation encountered in 
high-temperature superconductors in the experiments of p]. In such a case, 
the fermions are electrically charged excitations, which can represent excita- 
tions about the nodes of the d-wave superconducting gaps. The relativistic 
nature of the excitations is due to the nodal structure. These excitations may 
be holons, capable though of interlayer hopping, or could be real (nodal) elec- 
tron degrees of freedom. In this case the fluctuating gauge fields are assumed 
electromagnetic in origin. 

The fully relativistic case, examined in section 2, may admit a less con- 
ventional physical interpretation, and correspond to the case considered in 
[|I^. There, the abelian gauge interactions represent magnetic interactions 



in a spin-charge separation framework for high-temperature superconduc- 
tors, in which one admits intrasublattice hopping in the underlying doped 
antiferromagnetic model. In contrast to that work, however, which used 
non-relativistic fermions, here the fermions (holons) are relativistic, since 
one considers the spin-charge separation near the nodes. One also assumes 
a four-dimensional nature for such fermions, as a result of interplanar hop- 
ping. However, the above interpretation of the relativistic case is still not 
fully understood in the sense of not having been derived by an appropriate 
microscopic model at a satisfactory level of mathematical rigor. 

At any rate, the generic results obtained from the field-theoretic analy- 
sis of the present article may be useful if one wishes to compare them with 
various scenaria that may be in operation in realistic condensed-matter sit- 
uations, namely detailed phenomenological models derived from microscopic 
condensed-matter systems with relevance to doped antiferromagnets (and 
thus high-temperature superconductivity). In such cases it would be also 
useful to examine models where one has non-Abelian {SU{2) type) four- 
dimensional fluctuating gauge fields, coupled to fermions, in the presence of 
external magnetic fields. This would be more relevant to particle-hole sym- 
metric spin-charge separating doped models, such as those considered in [|]. 
We hope to come to a detailed discussion of such issues in a forthcoming 
publication. 
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A Anisotropic Schwinger-Dyson equation 

In this Appendix we discuss the derivation of the Schwinger-Dyson equation 
for the fermion propagator in the anisotropic case. To this end, we shall 



follow the derivation given in |[T^ for the isotropic case. Let us first recall 



some definitions. Starting from the Lagrangian 

C = -\t^,T^'' + ^ [z .^ - (7 4 - a; (z^s - ^^3) 7' - m] ^, (41) 
we define the connected graphs generator functional W by 



exp W% r], j^] = Jv[A^, ^, ^] exp {^ ^/ + ^ /^'''A + ^^ + ^^) 

' (42) 

W has the following functional derivatives (we do not write the space-time 
indices) 



6W 


- (U^) = lA^ 


sw 

5r} 


- (t^) = lip 


5r] 


- (1^) = li^ 


5rj 67] 


-ijip H /*'\1> 



(43) 



where the expectation value {O) of an operator O is 



(O) = lv[A,^,^ O exp h J C + i I [jA + T]^ + ^r])\ , (44) 
and we define 

Inverting the relations between {j^,f],'r]) and {A^,iIj,iIj), we define the effec- 
tive action T[A^,iIj,iIj] as the Legendre transform of W[jfj^,fi,ri] by 

W = iT + iJ (^fA^ + r]^ + ^ri) . (46) 

From this definition we extract the following functional derivatives: 



— = -f 



5A^ 



5r _ _ 

r ^ 

^# = -^ 

ir^ - -^--^(Iw^y. (47) 

^-^ 0-?/' (5-?/' yoT] or] J 

The starting point to derive the Schwinger- Dyson equation is to assume that 
the integral of a derivative vanishes, such that we can write 

lv[A^,^,^]-=—expl^tJc + tlfA, + ^r] + V'^] = 0, (48) 

which leads to, after a functional derivative with respect to rj and setting 
jt^ = 1] = f] = 0, 

S{zi - Z2) = {i (^^1 - m - ixd-i'f] ._. s ^ . 1 x lj>=»?=»?=o (49) 
^ ' dr]{zi} dr]{z2) 

■ ( ^ a ^ 3\ ^ TT. ^ , 

V?(5j^(2;i) «(5j^(2;i) / Sr]{zi) 6t]{z2) 
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We now have to turn this last equation into an equation for the effective 
action if we wish to obtain a relation between the proper functions. For this, 
we will use the last equation of Eq. ( ^7\ ) and we remark that 

■^r ^ f [ ^r SAM 6i^{y) ST T t Sjjjy) \ 

5j,{z) J y\5AM5j,{z) 5j,iz)6^{y) Sij{y) 5j,{z) j 



dy 



6r ( 6'^Y \ ^ 



bAM \bA^{z)bAM) 



\bA^{z)bi){y)) 6^iy) 



(50) 



6^{y) \6A,{z)S^{y)) 
such that we have for vanishing sources 

^ 5^{z{)5A,{y)5^{z^) \6A^{z2)6AMI ' ^ ^^ 

The proper functions (respectively fermion propagator, photon propagator 
and vertex) are defined by 

6^T 

Dj{Zi,Z2) 



,u v-1, -i; 5A^^{z^)5A-{z2) '^=^^=^-=0 

1 ^^r ^ 



k^{Zi]Z2,Z^) 



g 6A^{zi)6i){z2)6i){z:i] 



IV=5/'=Am=o 



such that we finally obtain from Eg . (^91) after multiplying by S (we do not 
write the space-time indices) 

G = S-A7ra f SYGA^GD^^ + xAna f S-f'^GA^GD^^ (53) 

where g^ = Ana. 
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